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ABSTRACT 

The study of Wolf-Rayet stars plays an important role in evolutionary theories 
of massive stars. Among these objects, ~ 20% are known to be in binary systems 
and can therefore be used for the mass determination of these stars. Most of these 
systems are not spatially resolved and spectral lines can be used to constrain the 
orbital parameters. However, part of the emission may originate in the interaction 
zone between the stellar winds, modifying the line profiles and thus challenging us to 
use different models to interpret them. In this work, we analyzed the HenA4686A + 
CivA4658A blended lines of WR30a (W04+05) assuming that part of the emission 
originate in the wind-wind interaction zone. In fact, this line presents a quiescent base 
profile, attributed to the WO wind, and a superposed excess, which varies with the 
orbital phase along the 4.6 day period. Under these assumptions, we were able to fit 
the excess spectral line profile and central velocity for all phases, except for the longest 
wavelengths, where a spectral line with constant velocity seems to be present. The fit 
parameters provide the eccentricity and inclination of the binary orbit, from which it 
is possible to constrain the stellar masses. 
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1 INTRODUCTION 

Massive stars are known to drive strong winds, which are 
responsible for the transfer of a large amount of the stellar 
mass to the interstellar medium, contributing to the feed- 
back of chemical elements, and to the creation of cloud cav- 
ities in which these objects are found. Typically, O stars 
present mass-loss rates of Mq ~ 10 -6 — 10 -5 Mq yr _1 and 
wind velocities of vq ~ 2000 — 3500 km s _1 , while Wolf- 
Rayet (WR) stars present M WR ~ 10~ 6 - 10~ 4 M yr" 1 
and uwh ~ 1000 - 4000 km s _1 (Nugis & Lamers 2000, 
Lamers 2001). 

In massive binary systems in which both stars present 
high mass-loss rates and high-velocity winds, the collision of 
the winds will occur. Therefore, a contact surface is formed 
where the momenta of the two winds are equal, surrounded 
by two shocks. The post-shocked gas, cools as it flows along 
the contact surface and is responsible for strong free-free 
emission at X-ray and radio wavelengths. X-rays in massive 
binary systems are orders of magnitude higher than those 
observed in single massive stars, and are used as an indica- 
tion of binarity. 

UV and optical lines can also indicate the binary nature 
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of a given object, when they present periodic profile varia- 
tions. If the lines are of photospheric or atmospheric ori- 
gin, as the stars move along their orbit, they suffer Doppler 
shifts, which depend on the orbital phase and inclination. 
However, some massive objects show periodic variable line 
profiles that cannot be explained under these assumptions. 

Seggewiss (1974) noted that the binary system WR79 
presented two peaks, superimposed to the Cm emission line, 
which changed their position and intensity with time. Typi- 
cally, in double line spectroscopic binaries, each peak moves 
in a different direction, indicating opposite velocity compo- 
nents along the line of sight for each star; however in WR79 
both peaks moved in the same direction. Liihrs (1997) pre- 
sented a model in which the two peaks were not produced by 
the stellar photosphere, but were generated by the flowing 
gas at the contact surface between the two strong shocks. 
This model reproduced well the data for WR 79, but failed 
to reproduce the line profiles of other WR binary systems, 
among them WR 30a (Bartzakos, Moffat & Niemela 2001). 
Falceta-Gongalves, Abraham & Jatenco-Pereira (2006) im- 
proved Liihrs' model introducing more realistic parameters, 
as stream turbulence and gas opacity, to account for the line 
broadening and peak displacement. 

In the present work, we applied this model to WR 30a, 
which was classified as a W04+05 binary system (Moffat 
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& Seggewiss 1984, Crowther et al. 1998); its binary nature 
was reported by Niemela (1995) based on spectral-line radial 
velocities obtained with a high temporal resolution. Later, 
Gosset et al. (2001) presented a detailed analysis of the spec- 
tra of WR30a for several epochs and were able to confirm the 
binary hypothesis and to determine the period of P ~ 4.6 
days. They also noted strong line-profile variations, which 
made more difficult the determination of the stellar mass ra- 
tio and the orbital inclination from standard methods. They 
concluded that the CIVA4658A line-profile variations were 
related to wind- wind collision processes, but did not model 
the lines under such an assumption. The same conclusion 
was reached by Bartzakos, Moffat & Niemela (2001) using 
the ClvA580lA line. Also, Paardekooper et al. (2003) pre- 
sented photometric measurements at V and B bands; the 
light curves confirmed the period obtained by Gosset et al. 
(2001), but also showed higher frequency variability in the 
V band, with timescale of hours. They concluded that this 
could be due to the strong variability of the CrvA580lA line, 
possibly related to the wind-wind interaction. 

In the present work, we tested the wind-wind shock 
emission hypothesis on the CIVA4658 A excess line pro- 
file variations measured by Gosset et al. (2001) using 
the model developed by Falceta-Gongalves, Abraham & 
Jatenco-Pereira (2006), which is briefly described in Section 
2. In Section 3, we show the results obtained for WR 30a 
and present a brief discussion, followed by the conclusions 
in Section 4. 



2 WIND- WIND EMISSION MODEL 

In the proposed situation, in which both stars present high 
mass-loss rates in supersonic winds, the contact surface, 
which is schematically shown in Figure 1, will have a ge- 
ometry described analytically by (Luo, McCray & Mac-Low 
1990): 



dy^ (r ] - 1 / 2 d 2 2 +d 1 2 )y 
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where D is the distance between the stars; d\ and d 2 are the 
distances of the primary and secondary stars to the contact 
surface, respectively, and r\ — M s v s /M p v p , where M p and 
M s are the mass-loss rates of the primary and the secondary 
stars, and v p and v a their respective wind velocities. The 
contact surface will asymptotically have a conical shape with 
an opening angle defined by /3, given by: 
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and the apex will occur at a distance to the primary star 
given by: 

dl|apcx = J D/(l+ J? 1/2 ). 



(3) 



Two shock fronts will be formed on both sides of the 
contact surface, generated by each wind, and the gas in the 
post-shock region will flow away along the contact surface. 
While it flows, the gas will cool due to expansion and radi- 
ation. This emission, detectable from radio wavelengths to 
X-rays, is the signature of wind-wind collisions (Usov 1992, 
Falceta-Gongalves, Jatenco-Pereira & Abraham 2005, Abra- 
ham et al. 2005, Pittard & Dougherty 2006). The stream of 




Figure 1. Schematic view of the wind-wind interaction surface 
and the model geometry. 



hot gas flowing along the contact surface will, eventually, 
reach temperatures that allow the recombination of certain 
elements. In this situation, the observer will detect the emis- 
sion line shifted due to the stream velocity component along 
the line of sight, which may intercept regions with differ- 
ent velocity components, as shown in Figure 1. The shaded 
regions in the shock layer represent the emitting fluid ele- 
ments. Here, observer A would detect two emission lines, a 
blueshifted line from fluid element 1, and a redshifted line 
from element 2. On the other hand, observer B would detect 
two blueshifted lines. 

Liihrs (1997) used these concepts to model line-profile 
variability during the orbital period of WR79. To reproduce 
the observations, he assumed a large emission region, lim- 
ited by two cones with aperture angles ft and ft + Aft, se- 
lected ad hoc to fit the the observations. However, the as- 
sumption is valid only in the case of quasi-adiabatic shocks, 
which occur in long-period systems. Radiative shocks evolve 
much faster, creating layers that are narrow and turbu- 
lent (Stevens, Blondin & Pollock 1992), and therefore can- 
not simply be described by a fixed range of beta values as 
assumed by Liihrs (1997). Falceta-Gongalves, Abraham & 
Jatenco-Pereira (2006) by-passed this problem by including 
turbulence in the stream layer, and also took into account 
the shocked gas opacity, which may result in line profile 
asymmetries. In that work, the line profile is obtained by 
integrating the emission intensities of each fluid element as: 
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where: 



fobs = v flow ( — cos ft cos if sin i + sin ft cos a sin ip sin i 

— sin ft sin a cos i), (5) 

is the observed stream velocity component projected into 
the line of sight, a is the azimuthal angle of the shock cone, 
i is the orbital inclination, ip is the orbital phase angle, a 
is the turbulent velocity amplitude, r is the optical depth 
along the line of sight and C{ip) is a normalization constant. 

In Figure 2, we show the model results for different or- 
bital and stream parameters. The dashed line model (a) was 
calculated using ft = 55°, i = 60°, a — 0.1 and Wfl ow = 1000 
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km s _1 ; the solid line (b) using (3 — 55° , i = 60°, a — 0.5 and 
»jo„ = 1500 km s" 1 ; and the dotted line (c) using (3 = 70°, 
i = 60°, a = 0.3 and «fl ow = 1000 km s _1 . It is noticeable 
that, e.g. in models (a) and (c), there are two peaks gen- 
erated by each of the shock-cone layers intercepted by the 
line of sight. As shown by Falcet a- Gonial ves et al. (2006), 
this is valid for 0° ^ i < 90°, while for i = 90°, the pro- 
file will present a single peak at tp = 0°, which corresponds 
to the cone axis coincident with the line of sight. Another 
comparative analysis between (a) and (c) shows the effect 
of changing (3; as f3 increases, the distance between the two 
peaks becomes larger. Models (b) and (c) present high tur- 
bulence amplitude, which leads to a larger line broadening. 
Here, we neglected the gas opacity (i.e. r << 1) to show the 
presence of two peaks in the line profiles. If the shock cone 
is optically thick, part of the emission along the line of sight 
will be absorbed, and the redder peak will be less intense 
or, eventually, completely disappear. 

In the next section we present the application of this 
model to the CIVA4658A excess line-profiles of WR30a, in 
order to determine the orbital and stellar parameters. 



3 THE CASE OF WR30A 

Gosset et al. (2001) presented a detailed spectroscopic study 
of WR30a with data covering 30 days, which corresponds to 
about 6 complete orbital cycles. From their observations, the 



strongest broad emission line is identified as a blend dom- 
inated by CIVA4658A and HeIlA4686A. Subtracting a con- 
stant parabolic base line, attributed to the WR stellar wind 
emission, they obtained the averaged residual spectra shown 
in Figure 3. This excess emission presents anomalous profiles 
and variability in both line intensity and central wavelenght, 
all related to the orbital phase. This effect could originate in 
selective wind eclipses, resulting in different profiles as the 
O-star is located behind or in front of the WR star along the 
orbital motion. However, this is not the case in WR30a be- 
cause due to its low orbital inclination no phase-dependent 
atmospheric absorption would be expected. 



3.1 The orbital parameters 

As it is clear from Figure 2, a wind- wind interaction model is 
compatible with the line profiles observed in WR30a. Hence, 
we applied the model described in the previous section to de- 
termine the binary-system parameters. An important model 
parameter in equation (4) is the optical depth r(a) of the 
absorbing material; Falceta-Gongalves et al. (2006) assumed 
that it is the absorption by the dense material that accumu- 
lates along the contact surface that produces the asymme- 
tries in the line profiles, while any other absorption is taken 
into account by the multiplying factor C(ip). Since most of 
the lines seem to be symmetric, at least within the uncer- 
tainties due to the subtraction of a phase independent base 
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Figure 4. Best fitting synthetic line profiles obtained for (3 = 50°, 
i = 20°, a = 0.3 and v fiow = 2200km s _1 . The orbital phase of 
each data set is exactly the same shown in Figure 2. 

profile, we assumed r(a) = 0. It turned out from our fitting 
that C(ip) is also independent of the orbital phase.. 

In order to model the observed excess emission spectra 
for each orbital phase, we performed calculations varying the 
opening angle, the turbulence amplitude and the stream ve- 
locity for different values of the orbital inclination. We were 
not able to fit simultaneously the complete width of the 
observed profiles and their mean velocities for all epochs, 
even for high values of the orbital inclination. Instead, wc 
obtained a better agreement when we considered that the 
emission at the longest wavelengths, seen as a peak near 
4690A is not formed in the shock region. The best fit oc- 
curred for (3 = 50° ± 5°, i = 20° ± 5°, v aow = 2200 ± 500 
km s~\ and a = 0.3 ± 0.1; the corresponding profiles are 
shown in Figure 4. The phase angles used for the fit for each 
epoch, and their respective uncertainties are shown in Figure 
5, as a function of the orbital phase. The uncertainties in the 
parameters were obtained by changing each one of them in- 
dependently until the model became incompatible with the 
observations. The excess profile for orbital phase 0.7 is very 
well-fitted by the wind-wind emission model, with the excep- 
tion of the already mentioned additional peak near 4690A. 
For orbital phases near 1.0 the model profiles are wider than 
those observed, which seemed to be single peaked. This may 
be related to the subtracted base line profile, which was ob- 



tained by Gosset et al. (2001) for phase 0.55 and equally 
applied to all orbital phases. 

If the subtracted profile is produced in the WR stellar 
wind, as it seems to be the case, it would be necessary to 
take into account the fact that the wind is not completely 
symmetric but shaped by the conic hole with aperture f3, 
as the wind is deflected at the shock. Since the orientation 
of this hole changes with orbital phase, we would expect 
also changes in the broad profile produced by the WR wind. 
For the particular case of i = 90°, when the WR star is in 
conjunction regarding the observer, there will be a stronger 
deficit of emission at the redshifted part of the line pro- 
file. On the other hand, if the WR star is in opposition, 
the amount of upcoming gas will be lower, affecting more 
the blueshifted emission. For i — 0° , the profile would not 
change as the orientation of the cone is always the same with 
respect to the observer. For WR30a, with i ~ 20°, at WR 
conjunction there will be mainly an effect on the red side 
of the profile but the blues side would also be affected. The 
opposite would occur at WR opposition. This is compatible 
with the narrower excess profiles near phase 1.0 as can be 
seen in Figure 4. 

At long wavelengths, the excess near 4690A could pos- 
sibly be the blended emission from the HeIlA4686A line; 
although it seems decoupled from the shock emission during 
the orbital period, it is too narrow ~ 20A (1300 km s _1 ) to 
be attributed to the WR stellar wind (> 2000km s _1 ). It 
could be, however, a fraction of the WR P-Cygni profile or 
it could also be related to the remnant shocked material sur- 
rounding the stellar system, ionized by the stellar radiation. 
Gosset et al. (2001) also attributed the narrow absorption 
superposed on the broad excess emission to a HeiiA4686A 
absorption line from the O star. 

The phase angles derived from our model provide addi- 
tional information regarding the orbit, namely its eccentric- 
ity, since the model strictly depends on the angle between 
the cone axis and the line of the sight, while the observed 
profiles depend on the orbital phase. The dependence be- 
tween the two parameters is related to the orbital eccentric- 
ity, as shown in Figure 5 for e = 0.0, 0.2 and 0.5. The dotted 
line represent the result for e = 0.0 (circular orbit) , in which 
the phase angle changes linearly with time. The dashed line 
shows the fast phase angle variation during the periastron 
passage for e = 0.5. The data for WR30a seem to fit best the 
curve corresponding to e = 0.2 when the lines are shifted by 
ipo = —20° in phase angle (0.056 of the orbital phase). Part 
of this shift could by explained by the Coriolis Effect, which 
produces a deviation of the cone axis as a consequence of the 
relative motion of the stars; however, Gosset et. al. (2001) 
estimated this effect as only 0.01 in orbital phase. Most of 
the contribution is probably produced by the binning of the 
data in intervals of one tenth of the period, which could also 
explain why the eccentricity we found is different from zero. 

3.2 The stellar masses 

After modeling the wind-wind shock emission, and deter- 
mining the orbital parameters of the binary system, we can 
constrain the stellar masses. For this purpose, we analyzed 
the velocity curve of the O-star obtained by Gosset et al. 
(2001) from measurements of four lines (HeIlA4542A, H7, 
HeIlA4200A and H8) during the orbital period. Their re- 
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Figure 5. The dependence of the phase angle on the orbital phase 
for eccentricities e = 0.0 (dotted), 0.2 (solid) and 0.5 (dashed), 
where <po = —20°. The data for WR30a (squares) were obtained 
from the phase angles of the fits and the observed orbital phases 
of Figure 3. 



suits are shown in Figure 6, as well as the fit of the modeled 
curve for e = 0.2 and Ko = 25 km s _1 . As a result, the 
mass function: 
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gives f(m) = 0.0074 M . If we assume Mo = 40 M (60 
M ), using i — 20° obtained from the shock emission, we 
find Mwr = 7.5 M (9.7 M ), which gives a mass ratio 
Mwr/Mo ~ 0.2. Regarding the orbit, using the obtained 
mass function and orbital inclination, we find ao — 5.4 R 
and «wr — 30 R . 

These results are in agreement with those expected for 
the masses of O and WO stars 40 - 60 and 8 - 10 M Q , 
respectively (Howarth & Prinja 1989, de Marco & Schmutz 
1999); also, the distance between the stars is larger than the 
sum of the expected stellar radii: D > Ro + Rwr ~ 20 R 
(Schaerer, Schmutz & Grenon 1997). 

From Equation 2 and the value of the cone-opening an- 
gle (f3 = 50°) we calculated the wind momentum ratio and 
obtained r\ = 0.12, which is compatible with typical ob- 
served values in systems similar to WR30a (Lamers 2001). 
From Equation 3, we estimated the distance from the apex 
to the stars: for rj = 0.12 and the total separation D = 35.4 



25 R and d 



2 | apex 



10 R . Consider- 



R we obtained d\ \ 
ing a typical O-star, tfelapcx would be of the same order than 
the stellar radius, and the frontal part of its wind could be 
crushed onto the stellar surface. The shape of the contact 
surface would be spherical near the apex, affecting maybe 
the X-ray emission from the system (Usov 1992), but our 
model will still be valid, since the line emission originates at 
larger distances from the apex, on the contact surface. 



4 CONCLUSIONS 

In this work we presented an application of the wind-wind 
shock emission model proposed by Falceta-Gongalves et al. 
(2006) to the line profile variations of WR 30a. In this model, 
the wind-wind shock structure can be represented by a cone, 
along which the shocked material flows. This gas will emit 
radiation, cool down and eventually reach recombination 





50 




30 


(km/s 


10 






o 
o 


-10 


CD 




> 






-30 




-50 









7 \ / \ 






\ ° / 


\ 7^ 


NEb-o P 


o-o / ° 







0.5 

Orbital phase 



Figure 6. The O-star radial velocity as function of the orbital 
phase. The observed data (circles) were obtained from Gosset et 
al. (2001). The solid line represent the fit for Kq = 25 km s _1 
and e = 0.2. 



temperatures. The observed emission lines will then suffer 
Doppler shifts due to the stream velocity component along 
the line of sight. During the orbital movement, the cone po- 
sition will change, as well as the radial velocity, which will 
cause the line profile variations observed in several massive 
binary systems. 

Gosset et al. (2001) obtained detailed spectra of WR30a 
during more than 30 days. They determined the orbital pe- 
riod (P — 4.6d), and obtained the radial velocity curve for 
the O-star. Regarding the WR component, they found that 
the blended HeiiA4686A and CivA4658A lines showed a vari- 
able excess emission. In the present paper we modeled this 
variable emission, being able to reproduce the variations ex- 
cept for the red part of the profiles, which seemed to be 
unchanged in velocity and were probably generated in the 
stellar wind of the WR star instead of in the shock region. 

The best- fitting result was obtained for (5 = 50°, i = 
20°, a — 0.3 and Vfl ow = 2200km s _1 . Also, correlating the 
orbital phase with the modeled phase angle, it was possible 
to determine the orbital eccentricity as e = 0.2, similar to 
the value of 0.0 previously assumed ad hoc by other authors. 
Although both values lead to very small differences in the 
orbital shape, its value is important for the determination of 
the stellar masses and orbital separation between the stars. 
Using this eccentricity and orbital inclination to model the 
radial velocity curve for the O-star, we found Ko = 25 km 
s" 1 . If we assume Mo = 40 - 60 M , we find M W r = 
7.5 — 9.7 M Q , and the orbital major semi-axis ao — 5.4 R 
and awn — 30 R . 

The model showed itself a powerful tool for constraining 
the wind and orbital parameters of massive binary systems. 
The fits did not matched the data exactly at all epochs, but 
considering the difficulties of subtracting the emission excess 
from the original spectra, the general shape and peak posi- 
tion variations were well reproduced. We must state that 
this is a very simple approximation, taking into account the 
complexities found in such systems. Numerical simulations 
could give a more detailed analysis, and probably more ac- 
curate values for the model parameters in future works. 
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